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Abstract 



in 



An investigation of the spatial fluctuations and their manifestations 
in the vicinity of the quantum critical point within the framework of 
the renormalized 4> 4 theory is proposed. Relevant features are reported 
through the Ginzburg-Landau- Wilson (GLW)-based calculations, com- 
bined with an efficient non perturbative technique. Both the dimension 
and size, but also microscopic details of the system, leading to critical 
^ . behavior, and strongly deviating from the classical mean-field approach 

far from the thermodynamic limit, are taken into account. Further, the 
important role that harmonic and anharmonic fluctuations and finite- 
size effects can play in the determination of the characteristic properties 
of corresponding various systems, involving phase transitions and critical 
phenomena, is then discussed in detail with emphasis on the qualitative 
validity of the analysis. 
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and anharmonic fluctuations , Gaussian and non-Gaussian approximation, Hartree- 
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1 Introduction 



Nowadays the theory of phase transitions and critical phenomena seems to be 
well developed in general. It enables us to obtain both universal and non-universal 
properties for many model systems. However, there remains a number of unsolved 
important questions among which one can cite, for instance, the role of fluctuations 
on the stabilization of the physical and electronic systems, the fluctuations that 
induced or prevent phase transitions, etc. There also exist some physical systems 
having as a common characteristics the fact that the complex microscopic behavior 
underlies macroscopic effects. 

In simple cases the microscopic fluctuations average out when larger scales 
are considered, and the averaged quantities satisfy classical continuum equation. 
Hydrodynamics is a standard example of this, where atomic fluctuations average 
out and the classical hydrodynamic equations emerge. Unfortunately, there is a 
much more difficult class of problems where fluctuations persist out to macroscopic 
wavelengths, and fluctuations on all intermediate length scales are important too. 

In the last category are the problems of critical phenomena. The critical phe- 
nomena in the thermodynamic limit is characterized by the divergence of the 
correlation length £ near critical point. Nowadays, the experimental techniques 
have become so advanced that the correlation length £ can be pushed up to several 
thousand A and the samples under study become comparable with £. As a conse- 
quence, the effects of finite size of the samples on the critical phenomena become 
increasingly important. Generally such effects depend on the shape of the sample, 
the boundary condition, the dimension of the system and the number of compo- 
nents of the order parameter. In fully finite or quasi-one- dimensional systems the 
phase transition is smeared out, whereas in thin films of thickness L the critical 
temperature T C (L) is shifted with respect to the bulk T c . 

On the other hand, during the first half of the last century after the discovery 
of superconductivity the problem of fluctuation smearing of the superconducting 
transition was not even considered. In bulk samples of traditional superconductors 
the critical temperature T c sharply divides the superconducting and the normal 
phases. It is worth mentioning that such behavior of the physical characteris- 
tics of superconductors is in perfect agreement both with the Ginzburg-Landau 
(GL) phenomenological theory (1950) pQ and the BCS microscopic theory of su- 
perconductivity (1957) [2J. The characteristics of high temperature and organic 
superconductors, low dimensional and amorphous superconducting systems stud- 
ied today, strongly differ from those of the traditional superconductors discussed 
in textbooks. The transitions turn out to be much more smeared out. The ap- 
pearance of superconducting fluctuations above the critical temperature leads to 
precursor effects of the superconducting phase occurring while the system is still in 
the normal phase, sometimes far from T c . The conductivity, the heat capacity, the 
diamagnetic susceptibility, the sound attenuation, etc. may increase considerably 
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in the vicinity of the transition temperature [3J. 

Recently, we proposed a microscopic renormalized Gaussian approach to critical 
fluctuations in the GLW model and finite-size scaling to describe the phase behav- 
ior and critical phenomena in very varied systems [J]. Within this more rigorous 
approach the effects of fluctuations are examined beyond the standard Gaussian 
approach (SGA) [3 [6] in more rigorous detail, and we are able to establish the 
insufficiencies of the mean- field theory (MFT), and also estimate the width of the 
critical region where corrections to MFT are important. The approach allowed 
us to obtain the effective functional of the Gaussian GLW Hamiltonian (Hclw) 
expressed in terms of the collective variables without any procedure for further 
systematic improvement (e.g. by considering "higher order" terms). 

In this paper, we construct the H GLW (in the space of the fluctuating fields) 
devoted to the study of harmonic and anharmonic fluctuations and their mani- 
festations in the near-critical region for the GLW model. In the former case [I] 
we took into account only the Gaussian approximation and showed that the ob- 
tained results were better and in good agreement with the experience as those 
found within the framework of the SGA method. Here we generalize the approach 
taking into account the powers of field higher than the first one. We need to con- 
struct a renormalized GLW which includes the quartic term that also takes into 
account the microscopic details of the systems. The importance and role of the 4 
interaction were emphasized by Fisher and Wilson [7J [8j [9j [10] . Its presence allows 
a low-temperature behavior and its absence leads to the divergence of the Hamil- 
tonian. The main ingredient in the analysis is as follows: the Fourier transform 
of the 4> 6 field theory can be seen as an interaction between the Fourier compo- 
nents of the order parameter, which are in fact a combination of harmonic and 
anharmonic fluctuations modes of the order parameter. This makes it possible for 
the coefficients of the <fi 2 and 4 terms to be renormalized by the 4> 6 term coeffi- 
cient. The approach is a strategy for dealing with problems involving many length 
scales. The strategy is to tackle the problem in steps, one step for each length 
scale. In the case of critical phenomena, the problem is, technically, to carry out 
statistical averages over thermal fluctuations on all size scales. The method is to 
integrate out the fluctuations in sequence, starting with fluctuations on an atomic 
scale and then moving to successively larger scales until fluctuations on all scales 
have been averaged out. The integration of this deviation in the self-consistency 
on the entire spectrum of the lattice vibration, conducted to a very improved self 
consistent energy. This improvement of the self consistent problem by the pro- 
cessed harmonic and anharmonic fluctuations goes, for instance, to reflect itself 
on the thermodynamic quantities and electronic parameters, in the neighborhood 
of the phase transition, in particular the singularities as function of the control 
parameter, resulting notably in a substantial modification of the critical point and 
thermodynamic quantities. The resulting self consistent problem of this approach 
is solved analytically, what permits us to extract an effective theory, with notably 
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(V r 0(r)) 2 + F L 



a mean-field critical temperature renormalized by anharmonic fluctuations. The 
present work does not present any renormalization based method in the sense of a 
renormalization group theory approach but rather a self-consistent method (SCM) 
improving on Landau theory (LT). Apart from calculations of exponents and scal- 
ing functions, it is necessary to develop techniques for obtaining the corrections 
to the asymptotic critical behavior, in terms of a small number of non universal 
parameters which can be fit by experimental results on different materials. In 
this way, it is hoped that a more rigorous confrontation between experiments and 
theory can be achieved. 

2 Formulation 

The starting point of our investigation is the following Hqlw functional, 

d d r 

where is the Landau free energy given by 

F L = a(T)<ft 2 (r) + 6(T c )0 4 (r) + u^\r) + ■■■. (2) 

Here, (ft is the model order parameter characterizing the mode displacement. The 
quadratic coefficient a(T) = ao{T — T c ) vanishes linearly as the temperature ap- 
proaches the mean-field critical temperature T c . For stability purposes, the co- 
efficient constant Mo is chosen such that uq > 0. The LT of second order (fer- 
romagnetic) phase transitions, for example, amounts to postulating the existence 

of a development of type b(T c ) = close to T c , what is not right for the 

T=T C 

superconductors where the microscopic relations between the GL parameters are 
such as j^jhj = ^jy^) where ((x) is the Riemann zeta-function (C(3) = 1.202...), 
ks = 1.38 x 10~ 23 W s/K is Boltzmann's constant and v a parameter which will 
be defined later. But just like the quadratic coefficient, the constant b(T c ) must 
be to improve to take into account both the microscopic details, dimension and 
size of the system. The real expressions of a(T) and b(T) within the framework 
of a renormalized theory will be defined later. £ is the coherence length of the 
sample (interpolated down to T = K). In superconductors, average extension £ 
of a Cooper pair as a measure of the distance within which the correlation form- 
ing Cooper pairs is active. In magnetic systems, £o represents the lowest length 
(microscopic) for spin correlations (spin waves), while in solids, it represents that 
of the vibrational or phonon modes. It is in general of the order of the interac- 
tion range. However, fluctuations with wavelength £o will be seen to be always 
negligible, d is the dimensionality of physical space. 

The polynomial in Eq. (2) originates from a power series expansion of some 
potential V(<ft) [TTj IT2"]. The order to which its terms are kept depends on physical 
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considerations given by renormalization theory [13] and on the type of symmetry- 
breaking effect [H] Eq. (2) is supposed to model. The odd powers in Eq. (2) are 
dropped as a result of e.g. time-reversal invariance. A "0 6 " model is of interest 
since, on contrast to "0 4 " models it can describe both second (b(T) > 0) and first 
order (b(T) < 0) phase transitions. Moreover, it displays a "butterfly" catastrophe 
[T5] which is more complex than the cusp catastrophe of the "0 4 " model. The "0 6 " 
model allows the ground state to be up to triply degenerate; situations with no 
degeneracy, double degeneracy or triple degeneracy can be studied by varying 
parameters in the expansion. The case a(T) = b(T) = is of particular interest in 
the context of tricritical points on phase diagrams [TS1 [T7] . 

LT [TTJ implicitly assumes that analyticity is maintained as all space-dependent 
fluctuations are averaged out. The loss of analyticity arises only when averaging 
over the values of the overall average (order parameter) <p. It is this overall aver- 
aging, over exp(— (3Fl), which leads to the rule that Fl must be minimized over 
(p. At temperatures that are far below, or far above a critical point, the behavior 
of the order parameter reassembles a tranquil ocean with no significant amount 
of thermal noise in its fluctuations. But fluctuations become increasingly impor- 
tant near the critical point as the correlation length diverges. At the second-order 
phase transition, infinitely long-range "critical fluctuations" develop in the order 
parameter. The study of these fluctuations requires that we go beyond LT. 

The SGA to the problem posed by the Hqlw is to decompose 0(r) into its 
Fourier components 0(q) according to 



According to Eq. (3) the limit on wavelengths means that the integration over q is 
restricted to values of q with \q\ < A. Averaging over long- wavelength fluctuations 
now reduce to integrating over the variables <p(q), for all \q\ < A [T5]. There are 
many such variables; normally this would lead to many coupled integrals to carry 
out, a hopeless task. Considerable simplifications will be made below in order to 
carry out these integrations. 

We need an integrand for these integrations. The integrand is a constrained 
sum of the Boltzmann factor ks over all atomic configurations. The constraints 
are that all \q\ < A are held fixed, this is a generalization of the constrained sum 
in the LT. We shall assume Landau's analysis is still valid for the form of HqysN) 
that is, Hqyw is given by Eq. (1). However, the importance of long- wavelength 
fluctuations means that the parameters a(T) and b(T) depend on A and then 
on the dimension d of the system. The d dependence of a(T) and b(T) will be 
determined shortly. However the breakdown of analyticity at the critical point is 
a simple consequence of this d dependence. Details will also be discussed shortly. 

The change of variables to the Fourier modes, the -£/glw expansion is given by 




(3) 



| 5 |<A 
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the expression 



g<A 

E 0(q / )^(q")0(q w )0(-q / - q" - q" 



{q}<A 

+ j^d E ^(qi)^(q2)^(q3)0(q4)0(q 5 )0(-qi - q 2 - q 3 - q 4 - q 5 ) ( 4 ) 

{g*}<A 

+ ••• 

where the g-mode function (5 _1 (q) = a(T) + cq 2 , and L is the linear dimension 
of the sample. Therefore, it is useful to express the partition function Z as a 
functional integral of the wave vector fluctuations 0(q). Accordingly, Z generalizes 
and therefore factorizes into 



Z = l D<p ■ exp 



+ ]7 E ^WW'¥-q' - q" - q" 

M<A 



w 



E ^(qi)^(q2)0(q3)0(q4)0(q 5 )0(-qi - q 2 - q 3 - q 4 - q 5 ) + 



te}<A 

(5) 

where Dip is given by 

P0 = IJ(27r)- 1 d0(q)^(-q), (6) 

and j3 = 1/ksT. T><fi is used to denote the measure of the functional integral. 
This quantity ensures that the total probability is normalized to unity through a 
constant of proportionality 11 (2tt) [131 E3- i^ 71 ) 1 left ou ^ is formally divergent 
in the thermodynamic limit; it does not affect averages that are obtained from 
derivatives of such integrals. 

In papers [20] and [21], a yet another approach to critical fluctuations has 
been proposed. It is based on the fact that quartic term is dominant and near 
criticality since a(T) — > as T — > T c . In the first paper [20] using mean-field 
approximation (MFA), the influence of homogeneous fluctuations was examined 
through the expansion 

Hmfa = #glw(0) + A 2 (0 - W + A 4 (0 - 0) 4 + ■ ■ • (7) 

where 4> is the equilibrium value of the order parameter and -£/glw is given by Eq. 
(1) without the Ginzburg term. In the second paper [2T] expansion like that in 
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Eq. (7) was studied for case other than spontaneous second-order transitions, i.e., 
for field-induced transitions, first-order transitions, and liquid-vapor transitions. 
In both papers calculations were performed using non-Gaussian integral given 
by[22j[23] 

poo 

/ <p 2mp ~ l exp (-A 2 2m - A 4 4m ) d<f> 
Jo 

= (2m)- 1 (2A 4 )^/ 2 r(p) J D_ p (A 2 / v / 2A;) exp (a 2 /8A 4 ), (8) 

and analyzed in all regimes including finite-size. D_ p is the parabolic cylinder 
function [23] and T is usually so-called gamma function. The most interesting con- 
clusions were that Gaussian approximation fails for all values of parameters except 
for in the thermodynamic limit (V — > oo). As a result, power-law predictions of 
the Gaussian prediction are incorrect, except for at V = oo, and should therefore 
be replaced by exponential asymptotic behavior as predicted by non-Gaussian 
methods. This results, for example, in vastly different asymptotic predictions 
for finite-size scaling as was discussed at length in both papers [201 [21] . These 
two papers, however, dealt exclusively with mean-field properties of non-Gaussian 
critical fluctuations, setting the Ginzburg term to zero. This, of course, neglected 
a very important property of critical systems, i.e., their spatial inhomogeneity 
PP. Nonetheless, interesting results were obtained [201 EI], leading asymptotically 
to the Gaussian approximation results and also providing finite-size scaling for 
T 7^ T c . Therefore, we believe that the proposed non-Gaussian method offers large 
region of analyticity and possesses better convergence properties. The present pa- 
per is intended to provide another insight into the problem by using non- Gaussian 
way of averages calculation and keeping a significant part of the Ginzburg term. 
We believe the renormalized Hqlw is capable of properly describing both the crit- 
ical and near-critical regimes. A non-Gaussian method of calculation, however, 
must be employed to adequately reveal the deviation from asymptotic behavior 
and thus the crossover phenomena. 

2.1 Non-Gaussian Fluctuations 

A major limitation of the GLW theory is its incapacity to account for the 
characteristics of the intermediate mode between the adiabatic mode on the one 
hand and the non-adiabatic mode on the other hand. It is also and especially 
its incapacity to take account of the microscopic details and the dimensionality 
of the systems. For instance, the GLW theory also predicts incorrect results, 
like an unphysical phase transition in ID, or incorrect critical indices, in higher 
dimensions. However, Scalapino, Sears and Ferrell [25J showed that, at least in ID, 
this failure is due to an improper treatment of fluctuations, so the complete success 
of the GLW theory will depend essentially on how fluctuations are taking into 
account. Both mean-field theories and Gaussian-based power expansions as well 
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as renormalisation-group calculations are concerned with the asymptotic properties 
of critical systems in the sense of infinite size and immediate proximity to T c . Thus 
they predict singular power-law behavior of the systems with universal exponents 
and scaling functions. As has been recently made abundantly clear, the asymptotic 
regime is indeed very small and agreement with experimental data deteriorates 
very rapidly outside the near-critical region; in fact extrapolation from the critical- 
state equations fails outside criticality and vice-versa, and the analytic noncritical 
equations of state do not reproduce the correct singular behavior at criticality. 
Thus there is a real challenge and a need to develop a method of calculation that 
both incorporates the asymptotic critical behavior and the crossover to regular 
regime. 

In the approximation that we propose, we wish to include those mode-mode 
coupling terms which involve balanced pairs of q and -q wave vectors, assuming 
that the remaining combinations are less important as they may lead to numerous 
cancelations. We can analytically understand the effective free energy of the expo- 
nential argument in Eq. (5) as that containing the quadratic term for the modes 
of free fluctuations of the order parameter <fi(q) (quadratic form), the quartic term 
for the modes of harmonic fluctuations of the order parameter <f)(q), and the sixth- 
order term as the anharmonic interactions between these modes. This condition 
seems to be effectively justified if the sixth-order coefficient of the expansion is 
small. 

Accordingly, it is important in this case to precise the following scenario: 
(i) The 2 and 4 terms are linked to the essential fluctuations, (ii) The <j) 6 
term is linked to the redundant fluctuations, (iii) All other terms correspond to 
the unessential fluctuations, i.e. contain no essential new physics and in fact are 
"irrelevant" (to the zero-temperature critical behavior in the sense described by 
Wilson PEIIH]). 

Therefore, the sixth-order term appears as and seems to play the role of inter- 
action terms between the Fourier components of the order parameter, which are in 
fact the combination of harmonic and anharmonic fluctuation modes of the order 
parameter. The sixth-order term can thus contribute to improve the quadratic and 
quartic coefficients. This makes it possible for the coefficients of the 2 term to 
be renormalized " an-harmonically" and the 4 term to be renormalized " harmon- 
ically" by the </> 6 term coefficient. It should keep in mind that the GLW theory is 
a phenomenological theory, founded on the intuition, with its own laws and rules. 
It accounts for the phenomena and it is that its justification. 

In this study, we will go beyond the GLW theory in order to evaluate the 
influence of the fluctuations on the critical phenomena and particularly critical 
temperature for systems with short-range interactions. We use a Hartree-Fock 
decoupling for 6 interactions in the consideration of the absence of long-range 
interactions. This approximation will consist in considering that the Fourier com- 
ponents interact (in harmonical and anharmonical manner) only through the mean 
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field produced by other modes. In the case of continuous phase transition this idea 
allows us to decouple the 6 term into a sum of product of two quantities with 
even powers: 

jr d E < ^ , ( c ii)^( c i2)0(q3)0(q4)0(q5)0(-qi - q 2 - q 3 - q 4 - q 5 ) 

{?i}<A 



* 15 E (jj* E (k(q / )0(q ,/ )0(q ,/, )0(-q / - q" - q"'))|) W(qM-<0 

q V {q t }<A J 

E ( j? E (k(q)0(-q)|) ) W) W)^(-q' - q" - q'")) + • <o) 

The factor 15 takes into account all possible contractions(statistical average of two 
or four modes among six). The terms between brackets act like mean fields. The 
sixth order term thus uncoupled consequently becomes a combination of quartic 
and quadratic terms, enabling us to obtain an effective non-Gaussian theory. With 
all these important preliminaries and considerations, one may then conveniently 
write the effective H GLW in the absence of an external field as 



d d r 



(V r 0(r)) 2 + a*(T)0 2 (r) + 6*(T))0 4 (r) + 



and the result of applying Eq. (3) to Eq. (10) is 



(10) 



<j<A 

E ^(q / )0(q // )0(q /// )0(-q / - q" - q'") + ■ ■ ■ (n) 

W<A 

where the renormalized g-mode function G* *(q) = a*(T) + cq 2 . 

The quadratic coefficient a*(T) and the quartic coefficient b*(T) are now re- 
spectively given by 

a*(T) = a(T) + Q d (T) (12) 
b*(T) = b(T c ) + Q d (T). (13) 

The term 

^d(T) = E (|0(q')0(q")<Kq"')0(-q' - q" - q'")|), (14) 

{%}<a 

includes a new contribution to the critical point and corresponds to the anharmonic 
variance of the order parameter at a single point in space evaluated at temper- 
ature T. This term competes with Landau quadratic coefficient a(T). Therefore 
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it determines the critical line and hence incorporates both the asymptotic critical 
behavior and the crossover to the regular regime. 

The renormalized quartic coefficient obtained within the renormalized GLW 
approach is now given by Eq. (13) with the correction term Qd(T). The term 



then also includes another new contribution to the quartic coefficient. This term 
competes with Landau quartic coefficient b(T c ) and defines the tricritical crossover 
exponent. Both scaling densities fld{T) and 0<f(T) are relevant, i.e., exhibit critical 
fluctuations. Now the quartic coefficient can be cancelled at some point highlight- 
ing the existence of a tricritical point. This fact represents the principal difference 
between the transitions corresponding to the tricritical Gaussian fixed point and 
the ordinary second-order Gaussian fixed point. Hence the tricritical fixed point 
can be characterized as the simultaneous instability of the system to two types of 
critical fluctuations. Correction to the molecular-field tricritical behavior due to 
critical fluctuations will be discussed in a forthcoming paper. It is found that the 
asymptotic tricritical form of some thermodynamic quantities or functions is not 
a power law but a power law multiplied by a fractional power of a logarithm. 

Taking a look at the influence of spatial fluctuations on the stability of the 
renormalized GLW analysis, the importance of fluctuations seems to be evident 
and the divergence of the critical coefficient of the second-order term a* (T) is now 
related to the system size L and dimensionality d of the system. Therefore it 
is useful to determine the temperature at the critical point where the divergence 
of quadratic coefficient should be observed. Accordingly, the modified critical 
temperature is given by 



where Qd c — &>d{T*) is the solution of the self-consistent Eq. (25) that will be 
discussed in the next subsection. 

Eq. (16) shows that T* can be neglected if the anharmonic fluctuations are 
too high. However, taking into consideration the 6 term contained in the -£/glw 
functional leads to the emergence of the limitations of the GLW approach and 
MFA related to the critical temperature T c , which is presented here as a charac- 
teristic scale of temperature related to both thermic fluctuations and finite-size 
effects rather than the transition temperature [3]. Although that limitation is a 
disadvantage to the exactness of the GLW approach or MFA, its predictions are 
not at all lacking interest and only the importance of the anharmonic term Q^c 
should determine its degree of validity. 

Structural phase transition (for example) is accompanied by a change in struc- 
ture. Some of these changes in structure occur without macroscopic diffusion of 
matter in solids. They are initiated by local motions of atoms or molecular groups 




(15) 



q 




(16) 
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which can distort the lattice in the high temperature phase to form structures 
of lower symmetry at lower temperatures. These movements around equilibrium 
positions do not occur instantaneously at temperature T c , and they are actually 
initiated at a temperature lower or higher than the transition temperature [26]. In 
a way, this pretransition or pretransformation phenomenon is the equivalent of the 
nucleation process. Eq. (16) can thus contribute to clarify this viewpoint, since 
Qdc which is the solution of a n degree polynomial can be a positive or negative 
quantity. 

2.2 Mean-field approximation of Cld(T), Qd(T) and thermodynamic quan- 



A correct treatment of f^(T) shows that @d{T) is much more complex. Once 
critical fluctuations are not treated as the constants, one could imagine expanding 
order parameter in a Taylor's series about its value at some central location r . 
This means that i^GLW could be a complicated functional of </>, an expression that 
is hard to write down, with several parameters, instead of the simple GLW form 
with only two parameters a(T) and b(T c ). Consequently, a natural recourse is to 
use Gaussian approximation. 

Gaussian measures play a central role in many fields: in probability theory as a 
consequence of the central limit theorem, in quantum mechanics, in quantum field 
theory, in the theory of phase transitions in statistical physics. In this section, 
calculations will be performed and analyzed in all regimes including finite sizes, 
using the Gaussian integrals in the form [23], 



We can formulate the Hamiltonian by averaging the new contributions. The pro- 
cedure for identifying the variational parameter of the quadratic and quartic coeffi- 
cients, fld{T) and Gd(T), and for determining the value of T* and thermodynamics 
quantities is as follows. First the Hamiltonian is derived in the MFA and the ob- 
tained O t j(T) is substituted for the self-consistent equation. The derived equation 
is used to determine both the value of T* and thermodynamic quantities. 

As a result, it is shown that the value of T* and thermodynamic quantities 
strongly deviate from the classical MFA but thermodynamic quantities adopt 
mean-field critical exponents provided by the classical Landau approach. We find 
that the variational parameter of the Hamiltonian obtained by renormalizing that 
of the quadratic term, corresponding to the square of angular frequency, is pro- 
portional to T* — T. 



tities 




p>0, n = 0,l,2,---; 



(2n - 1)!! = 1 • 3 • 5---(2n - 1). 



(17) 



n 



The thermal average of a physical quantity y(4>) is given by 



(y) 



y((p) exp[-(3H eS ]d(j) 




exp[-(3H eS 



(18) 



The obvious difficulty that the transformation given by Eq. (3) brings relates 
to the mode- mode coupling present in the last term in Eq. (11). Consequently, a 
natural recourse is to use the Gaussian approximation where H e g[(f)] is truncated 
to 

/J o [0] = ^G* _1 (q)0(q)0(-q) (19) 

q<A 

This then conveniently factorizes the partition function as 



2o=n 



7T 



9 <iy /3(a (T-T c ) + ft d + g 2 )' 
By making a transition to continuum through 

/oo 
(■■.)d d q, 

the free energy is 

k B TL d rco 



(20) 



(21) 



F 



2(2tt)' 



Jo 



In 



.p(ao(T-T c ) + n d + q< 



d d q 



(22) 



The heat capacity of the system is given by 

)'(■ 



where 



d 2 F 

or 2 



T 

± C J 



a T c J 



+ less singular terms, (23) 



x dx 
[1+x 2 ) 2 '- 

27T d / 2 



T -T r 



e = 



(24) 



2 d - l ir d / 2 a 2 T(d/2) J 

The element volume is $ldq d ~ 1 dq, with = r(rf/2) ^ ne s Ph ere um ^ surface in 
space. Kd is an integral correction constant which is lattice dimension dependent. 
The behavior of the integral correction constant changes dramatically at d — 4. 
For d > 4 the integral diverges at large x and is dominated by the upper cut- 
off A, while for d < 4, the integral is convergent in both limits. Although the 
expression of the specific heat is presented in not so complicated form, it is not 
quite transparent to know what type of behavior it will exhibit depending on 
dimension. Only the nature of fl d will make it possible to describe its behavior. 
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The dominant behavior of Cq close to T* is through Cq ~ |e*| rfi/ ~ 2 , where the 
reduced and renormalized temperature e* = e + and v = 1/2, so that we 
obtain a = 2 — dv as in the SGA . When thermal fluctuations can be neglected, 
as it appears to be the case in conventional superconductors, the specific heat 
exhibits at the transition temperature T C (L — > oo) a step discontinuity. This 
differs drastically from the behavior when thermal fluctuations dominate. Due 
to the finite-size effect, the specific heat peak occurs at a temperature T*(L) 
shifted from the homogeneous system by an amount proportional to L -1 ^, and 
the magnitude of the peak located at temperature T* (L) scales as L a l v [27] . This 
result seems to be qualitatively in good agreement with the well-known result since 
C ~ \T-T*\ du ~ 2 0. 



2.2.1 Expressions of fluctuating quantities il^T), Qd{T) 



The structural phase transition has been, hiterto, treated as follow [2"g]: In the 
" quasiharmonic" approximation [23 [3D], it was assumed that harmonic frequency 
for the unstable mode is purely imaginary, namely, Uq < under the short- and 
long-range force constant competing with each other. Then the contribution from 
the anharmonic term was derived self-consistently, which is proportional to T. As a 
result, the square of the soft-mode frequency was given by u 2 = u 2 + c{u 2 ) oc T — T c 
where u denote the displacement from the interatomic distance. Other soft-mode 
theories are LT and MFT [3T| 132] . In the expansion of free energy, the normal 
coordinate was adopted as the order parameter. It was shown that the critical 
coefficient of the second-order term, a(T), equals the square of the soft-mode 
frequency, a(T) = u 2 oc T c — T. The LT gives a qualitatively correct view of the 
soft mode, but cannot explain how it occurs; furthermore, this theory is not valid 
close to T c where critical fluctuations can no longer be neglected. So its complete 
success will depend essentially on how fluctuations are taken into account. 

In the present work, we show the approach to this problem from the micro- 
scopic point of view, both in terms of the order parameter and self-correlation 
functions (0 2 ) and (0 4 ). From Eqs. (17), (18) and (21), we obtain the result that 
the temperature dependence of the variational parameters, fld(T) and Qd(T), cor- 
responding to the anharmonic and harmonic fluctuations of the order parameter, 
are presented by 



Kd 
V 



T 



1 2 



e + 



a, 



aoT c J 



where V 



(25) 



and 



0,i = /c„ 



T 

Tr 



e + 



-i 



(26) 



aoT c J 

K d and /Q are the dimension- dependent constants, respectively, and explicitly 
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45 % ^y /- + °° x d - l dx 



established as follows 

= 2 d - 1 n d / 2 a 2 T(d/2) J (1 + x 
and 

= 15^^ /- +00 x d -Hx 
^ 2 d - 1 7r d / 2 a r(rf/2) i (1 + x 2 )' 1 ' 

The equations ((25) and (26)) are derived taking into account the expected fact 
that fid(T) and Qd{T) do not explicitly depend on wave number q. The behavior 
of Kd is the same one as that of k^. For certain values of d, the integral /Q diverges 
at large x. This is not a first time when such problem arises and we know how to 
deal with it: this ultra-violet (UV) divergence is related to the restrictions on the 
applicability of the GLW functional for \q\ > so the integral has to be cut off 
at £ • |g| = x c ~ 1. 

It is well established that all correlation functions (in dimension d < 4) have a 
large cut-off limit after a simple renormalization, that is, after one has taken the 
deviation from the critical temperature as a parameter. The field amplitude renor- 
malizations are finite. In the critical phenomena situation instead, the fluctuations, 
which is related to microscopic parameters of the theory, is fixed. This means that, 
after the introduction by rescaling of the cut-off A, quadratic coefficient remains 
finite when A — > oo for T ^ T*. 

As we are concerned with phase fluctuations, Eqs. (25) and (26) show that 
their effect is dimension d dependent. It has been demonstrated that the inclusion 
of phase fluctuations leads to a reduction in the degree of order in d > 2 and to 
its complete destruction on d < 2; for d > 2, the phase fluctuations are finite 
while they become asymptotically large for d < 2 [5l [6j [33]. The more general 
result has been established and known as the Mermin- Wagner theorem [31]: it 
states that there is no spontaneous breaking of a continuous symmetry in systems 
with short-range interaction in dimension d < 2, and as a corollary the borderline 
dimensionality of 2, known as the lower critical dimension di has to be treated 
carefully. With regards to Eq. (26) dimension d = 2 seems to play a crucial role in 
agreement with the Mermin- Wagner theorem. Indeed, when d = 2, the quantity 
d/2 — 1 vanishes and the self-consistent equation for the harmonic correction term 
Qd (Eq. (26)) is reduced to a linear function of T with a positive coefficient as 

©d = fcd w ■ On the other hand, above d = 4, the MFT predicts correctly the 

universal quantities, whereas it is definitely not valid for dimension 4 and below. 
(Eq. (25) seems to arise this aspect with the quantity d/2 — 2 which is canceled 
when d = 4, and the self-consistent equation for the anharmonic correction term 
fid (Eq. (25)) is reduced to a function of T 2 with a positive coefficient as Qd = 
r 1 2 

Kd y~ . Although the expressions of the correction terms are presented in not so 
complicated form, it is not quite transparent to know what type of behavior they 
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will exhibit depending on dimension, temperature and the size of the system. The 
possibility that 9^ or Q d increases with decreasing temperature could highlight 
the quantum character of the system according to both the dimensionality and 
the size. However, the numerical aspect which will be approached in another 
work, will enable us to better clarify the behavior of 0^ and Qd and to give them 
a suitable physical direction. 

The renormalized critical coefficient of the second-order term a*(T), equals the 
square of the renormalized soft-mode frequency, a*(T) = u* 2 oc T* - T, but the 
critical temperature now takes into account both the dimension of the system, its 
size and critical fluctuation. 



2.2.2 Thermodynamic quantities 



To determine the degree of validity of this approach, let us formulate the 
Ginzburg criterion [35, 36j, which usually tells us quantitatively when MFT is 
valid. It is clear that the fluctuations become more and more pronounced as the 
temperature approaches the true critical point T* . The Ginzburg criterion in- 
dicates in a semi- quantitative manner the temperature range where the distance 
from the SCM is important. 

Therefore, the critical Ginzburg width [35| I3~T] in the vicinity of the critical 
temperature is then given by 



At 



G 



dc 



aoT c 



(29) 



In the last formula, Q dc — ^d{T*) is the solution of the following self-consistent 
equation (obtained in the same spirit as Eq. (19) in |Q): 



d --2 
2 



In 



e + 



aoT c 



In 



K d (a T c 



fide) 2 



(30) 



Here, T* also determines the temperature under which the description of fluctua- 
tions goes beyond the independence of the MFA modes due to precise agreement. 
Ate corresponds to the width of the critical region about the real transition tem- 
perature for which both the Landau and the GL theories are valid. Further by 
taking into account the effective GLW Hamiltonian, this quantity also defines the 
domain of validity of the classical critical behavior [38]. The classical description 
fails for |At| = \T/T* — 1| <C |Atc|. This finding is corroborated by the Ginzburg 
criterion and can be interpreted as [131 EH] • Thus, the Ginzburg criterion allows 
us to restore some credibility to the MFT in those cases. Eq. (25) shows that the 
anharmonic fluctuation Q d — > when V — > oo involving a*(T) — > a(T). In the 
thermodynamic limit, the renormalized quadratic coefficient a*(T) is the same as 
in the LT. It is well established that power-law predictions of the LT are correct 
in the thermodynamic limit, T* — > T c when V — > oo. 
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As it is well established (for V ^ oo), no Gaussian approximation even renor- 
malized in some way can describe precisely critical phenomena near the transition 
point. Also, within this approach, the effects of spatial fluctuations are strongly 
dependent on lattice dimensionality and we can now appreciate how fluctuations 
and correlations modify the macroscopic thermodynamic properties. 

From the thermodynamic definition, the inverse renormalized susceptibility is 
given by the following analytic expression: 



X 



a (T - T c ) + a 



(31) 



Taking a look at the influence of spatial fluctuations on the stability of the renor- 
malized mean-field analysis, the importance of fluctuations seems to be evident 
and the divergence of the susceptibility is related to the system size L and dimen- 
sionality d. 

The correlation length gives information about the distance for which the order 
parameter varies in the space. In order to calculate the correlation range, we 
first evaluate d 2 G*(q)/dq 2 where G*(q) is the renormalized g-mode autocorrelation 
function. Taking into account Eq. (20), the function G*(q) is obtained as 



G*(q) = (|0(q)| 2 ) 



2P[a(T) + n i + <f 



(32) 



Accordingly, as T 



T* within the limit q 



0, it diverges according to 



G*(q,T = T*) ~ q v ~ 2 , with a small positive value rj and the critical exponent 
f] = 0. Therefore, G~ 1 (q) appears as the reduced inverse q-mode autocorrelation 
function. According to Eq. (32), the renormalized correlation length is obtained 

as 



Q(T) 



_ limG .-. (9) ffl 

q^o V ^ V dq 2 



1/2 



-1/2 



, for T > T* 



-1/2 



(33) 



for T < T* 



which generally sets the characteristic length scale of fluctuations. The universal 
exponents and amplitude ratios are again recovered from this equation with regard 
to the critical temperature T*. Hence, the corresponding critical exponents are 
v — v' — 1/2. The temperature dependence of the renormalized correlation length 
Q(T) based on the above formula depends on the dimensionality d and the size L 
of the physical space and behaves differently in function of them. 

For finite-size systems (in d = 2,3) it has been recognized 0139] that the system 
size L " scales" with the correlation length of the bulk system. Thus it is convenient 
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to define a reduced length of the system: 



I* = L/C = kT c 



r 1/2 fe-r + -(n dc -n d )l 

L a n -I 



for T < T*. 



(34) 



The value of /n ~ £ 1 depends only on the nature of the substance under study. In 



fact both I* and l are dimensionless, and they are the ratio of the real thickness 
of the system to certain characteristic length. Indeed, if I* >> 1, no significant 
finite-size effects should be observed. On the other hand, for I* < 1, the system 
size will cut-off long-distance correlations so that an appreciable finite-size round- 
ing of critical-point singularities is to be expected. This result is not surprising 
considering the physical meaning of the correlation length which can be regarded 
as an indication of the influence range of the boundary condition. One should 
notice that what plays a role is not the real length L of the system but the reduce 
one which depends also on the deviation from the bulk critical point. Details 
will be discussed later in the section reserved for the applications. This shows 
that if the reduced length of the system is too short, either due to the small size 
of the system or due to its closeness to the bulk critical point, the influence of 
the boundary will strongly dominate. As I* becomes large, the renormalized Hq 
becomes independent of details of the system at the atomic level. This leads to an 
explanation of the universality of critical behavior for different kinds of systems at 
the atomic level. Liquid-gas transitions, magnetic transitions, alloys transitions, 
etc., all show the same critical exponents experimentally; theoretically this can be 
understood from the hypothesis that the same "fixed point" interaction describes 
all these systems. 

In order to illustrate the meaning of the correlation length, it is perhaps worth 
calculating the space-dependent correlation function G* d (X) defined by the Fourier 
transform of the q-mode autocorrelation function as 



which can also be considered here in the sense of the entanglement in the model 
system. The propagator G* d (X) should always go to zero at large distances, so that 
a measure of order in the system is the long-length behaviour of the propagator. 
If the propagator goes to zero then the system can only have short-range order. If 
the propagator goes to a non-zero constant at large distances, then we must have 
a non-zero order parameter, and the system has long-range order. 

The integral (35) is well behaved at small q. At large q without the exponential 
the integral would only be convergent for d — 1. For d = 2 there would be a log 
divergence, while for higher d there would be an algebraic divergence. Clearly this 
divergence is cut off by the exponential factor. As the exponential goes to 1 when 
X goes to zero, the correlation function should "diverge" at small distances. In 
fact the divergence will be cut off by the discretization length A (discreteness of 




(35) 
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underline lattice). The choice of scale depends on the situation we want to study. 
First of all, in a numerical attempt to describe experimental data, it is clear that 
for most practical cases the cutoff A must not be chosen too large, since this would 
require a knowledge of the lattice dispersion way beyond the parabolic approxima- 
tion corresponding to the simple gradient correction in real space. However, for 
the study of phase transitions, this is not a severe restriction since the anomalies 
appearing in such transitions result from long range correlations, i.e. the behavior 
of the system at small g-vectors. For large distances X > £ , the exponential 
is oscillating much faster than any other variations, and the correlation function 
should fall rapidly to zero (in fact exponentially). This behaviour is demonstrated 
if we explicitly calculate the integral for dimensions one, two and three. 

Using Eq. (35) and due to the fact that the space-dependent correlation func- 
tion assumes various forms depending on the dimensionality of the physical space 
0, G* d (X) satisfies 



G* d (X) 



exp 



5F*o(£), 
k (l^) ex P 



for d 



for d 



for d 



(36) 



where p = (x 2 + y 2 ) 1 ^ 2 , and r = (x 2 + y 2 + z 2 ) 1 ^ 2 . Ko(r) is the modified Bessel 
function (that is logarithmic for small arguments and exponentially decaying for 
large arguments). 

To keep the discussion self-contained, in the remainder of this subsection we 
review the predictions for the static critical exponents. First, we define the reduced 
temperature e* = (T — T*)/T c . The exponents a, (3, 7, 77 and v describe the 
singular behavior of the theory with strictly zero renormalized quadratic coefficient 
as e* — y 0. For the specific heat, taking into account Eq. (23) one finds 



C(T) ~ I e*| a + less singular terms 
rj and v describe the behavior of the correlation length £, where 

G^(r)E(^)^(0y-(i)(^) 



(37) 



(3? 



and the exponent 77 is defined through the behavior of the Fourier transform of the 
correlation function: 

G aP {q^U)~q- 2+ \ (39) 
The correlation length exponent v is defined by 



/*- 1 4 
4 ~ e 



(40) 
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The exponent (3 will be denned later. The last exponent, 5, is related to the 
behavior of the system in a small magnetic field h which explicitly breaks the 
0(A) symmetry. The six critical exponents defined above are related by four 
scaling relations [6]. 

a = 2-du, a + (3(1 + 5) = 2 = a + 2/3 + 7, 7 = (3(5-1). (41) 

3 Renormalized 4 theory 

All perturbative approaches are based on the division of the free energy into a 
Gaussian term and higher order perturbative terms. In fact, the coupling constant 
with <j) 4 model of the perturbation is not necessary small, so that the convergence 
of the perturbation expansion cannot be ensured. Thus some more effective ap- 
proaches to the calculation are needed. For systems with boundaries, one should 
consider the influences of the boundaries on the thermal properties near the bulk 
critical point. In additions, the spatial distributions of the order parameter should 
be taken into account for finite-size systems. As is well known for finite-size sys- 
tem, however, the spatial distribution cannot be considered as uniform any longer 
due to the influence of the boundary though the condition of minimum free energy 
would prefer a smooth distribution. 

3.1 Continuous phase transitions 

Different thermodynamic phases are characterized by certain macroscopic, usu- 
ally extensive state variables called order parameters; examples are the magneti- 
zation in ferromagnetic systems, polarization in ferroelectrics, and the macroscop- 
ically occupied ground-state wave function for superfluids and superconductors. 
We shall henceforth set our order parameter to vanish in the high temperature 
disordered phase, and to assume a finite value in the low-temperature ordered 
phase. Landau's basic construction of a general mean-field description for phase 
transitions relies on an expansion of the free energy (density) in terms of the or- 
der parameter, naturally constrained by the symmetries of the physical system 
under consideration [40J. For example, consider a scalar order parameter ip with 
discrete inversion or Z2 symmetry that in the ordered phase may take either of 
two degenerate values ip± = ±\ipo I [HI- We shall see that the following generic 
expansion (with real and renormalized coefficients) indeed describes a continuous 
or second-order phase transition: 

F* L = a*(T)^ 2 + b*(T)^ + ... -hip, (42) 

if the temperature-dependent parameter a* changes sign at T* . For simplicity, and 
again in the spirit of a regular Taylor expansion, we let a*(T) = a^(T — T*), where 
T* denotes the critical temperature from mean-field approach with renormalized 4 
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model. The free energy is of almost the same as the Landau functional, except for 
the presence of the fluctuations-dependent in the quadratic and quartic coefficients, 
which contain the essential information about the microscopic nature of the system, 
its size and dimension. Our functional therefore describes a set of interacting, 
weakly-GLW-damped excitations. 

Stability requires that b* > (the size of the sample is such that the renor- 
malized quartic term is a positive quantity); near the renormalized critical point 
we can simply evaluate b* at T*. Details will be discussed shortly. Note that the 
external field h, thermodynamically conjugate to the order parameter, explicitly 
breaks the assumed Z 2 symmetry ip — > —if). Minimizing the free energy with re- 
spect to -0 then yields the thermodynamic ground state. Thus, from dFl/dip = 
we immediately infer the equation of state 

h(T, i/j) = 2a*tp + 46* V 3 (43) 

and the minimization or stability condition reads < d 2 F*L/dip 2 = 2a* + 12b*ijj 2 . 
At T = T*, Eq. (43) reduces to the critical isotherm h(T*,ip) = 4b*ip 3 . For 
a*{T) > 0, the spontaneous order parameter at zero external field h — vanishes; 
for a*(T) < 0, one obtains ip± = ±\tpo\, where 

i) = vVl/2&* = f(T) ■ 0o (44) 



J o = v l a l/26 is the order parameter in LT. The behavior of the renormalized 
order parameter defines j3: 

(H) ~ for e* < (45) 

Taking into account Eqs. (12) and (13), the size dependence of ip enters in the 
amplitude prefactor f(T) defined by 




The SCM (taking into account both the dimension, the finite-size effects and the 
temperature dependence of f(T)) reveals a competition between three scales of 
energies that are in competition: the thermal energy a(T) oc /c^T versus the en- 
ergy resulting to anharmonic instabilities Qd on one hand, and the anharmonic 
instabilities versus the harmonic instabilities on the other hand. These competi- 
tions determine both the existence of the transition, but also the adiabatic regime 
of fluctuations of the "pre-transition" . 

As in many case, it is also possible to write the generalized GLW functional in 
a form which preserves certain transformations of ip(r) like either of the three. 

reversal ^(r) — > —ip(r) (47) 
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change of phase ^(r) — > e* 9 ?/>(r) (48) 

rotation ^(r) ->■ C/^(r) (49) 

(U is a rotation matrix). This holds for example for the Ising model, superfluid 
helium and the Heisenberg model respectively. The amount of ip is uniquely de- 
fined, but the sign, phase and direction of ip(r), are not defined. It depends on the 
history (preparation in an external field which removes the symmetry Eq. (48) or 
accidental fluctuations) of the system. 

Note the emergence of characteristic power laws in the thermodynamic functions 
that describe the properties near the renormalized critical point located at T = T c *, 
h — 0. Inserting Eq. (44) into the Landau free energy Eq. (42) one finds for T < T* 
and h — 



= ^ = ~w ( 5 °) 

Because of our assumptions about a* and b* , the renormalized free energy is pro- 
portional to (T* — T) 2 . This is characteristic of all second-order transitions, and 
consequently for the specific heat 

(51) 

whereas by construction F£(0) = and C£ =0 = in the disordered phase. Thus, 
Landau's renormalized MFT also predicts a critical point discontinuity 



AC" 



h=0 




T=T* 



(52) 



for the specific heat. Experimentally, one indeed observes singularities in thermo- 
dynamic observables and power laws at continuous phase transitions, but often 
with critical exponents that differ from the above mean-field predictions. Indeed, 
the divergence of the order parameter susceptibility indicates violent fluctuations, 
inconsistent with any mean-field description that entirely neglects such fluctua- 
tions and correlations. Let us mention that the jump of the heat capacity was 
obtained because of the system volume was taken to infinity first, and after this 
the reduced temperature e* was set equal to zero. 

Within the framework of this approach, taking into account the fact that 

= and b*(T*) = b(T c ). the anomalous part of the specific heat is 

T=T* 

^ c 

really given by 

AC* h - = L d T c -^-(l-^-). (53) 
2b(T c ) y a T c J 1 ; 

Then, it can be concluded that (taking into account new fluctuating quantities) if 
the transition takes place at low temperatures, the anomalous part of the specific 



dn d 

dT 
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heat could be negligible or important compared to Debey's specific heat stemming 
from the acoustic phonon distribution [3J. Our approach suggests the first case 
where this anomalous part could be too small to be detected. Eq. (53) is an 
indication of the decreased anomalous part of the specific heat as the anharmonic 
fluctuation of the system increases. 

A close look at the heat capacity transition in optimally doped YBCO samples 
shows that it starts several degrees above T c , and presents a rather sharp peak, 
with an increasing slope (dC*/dT) as T c is approached [32]. These are strong 
indications that thermodynamical fluctuations are playing an important role in 
the transition. We can use Eq. (29) to estimate the width of the critical region in 
a typical cuprate. 

A different deviation from mean field behavior is also seen in the transition of 
high T c cuprates, as shown in [32]. In addition to some broadening observed above 
T c the shape of the main transition is modified. Instead of a jump, it looks more 
as a narrow peak. Qualitatively, it reminds of the specific heat peak seen at the 
transition of superfluid Helium. 

An even more radically different form of heat capacity transition is observed 
in Bi 2 Sr2CaCu20s+s- There is no more heat capacity jump at the transition, but 
rather a cusp. The transition is better fitted by a Bose-Einstein condensation than 
by a BCS one [33] ■ More details will be discussed shortly. 

4 Examples of applications 

In order to test the theory of critical phenomena it is important to have accurate 
experiments on well characterized systems, very close to the critical point. Such 
experiments exist in magnetism, ferroelectric thin films and superconductors. 

4.1 Magnetism 

Magnetism is caused at the atomic level by unpaired electron with magnetic 
moments, and in a ferro magnet, a pair of nearby electron with moments aligned 
has a lower energy than if the moments are antialigned [IB]. The Curie point of 
a ferromagnet will be used as a specific example of a critical point. Below the 
Curie temperature T c , the ideal ferromagnet exhibits spontaneous magnetization 
(<f) 7^ 0) in the absence of an external field; the direction of the magnetization 
depends on the history of the magnet. Above the Curie temperature, there is no 
spontaneous magnetization. This ferromagnet ism is observed in certain metals like 
iron, nickel and cobalt. Just below the Curie temperature the mean field magnetic 
susceptibility is observed to behave as 

C 

Xm = T7f — 7Fn> 7 = 1 (54) 
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that is the Curie- Weiss law. C is the Curie constant [HI H2]. Experimentally, 
we observe that 7 is about 4/3. On the other hand, the reduced magnetization 
(f) = cj)/(Nn B ) (that is also the local spin density) is observed to behave as 

x { (T c - T)/ 3 T <T C (55) 

i.e. the exponent /3 is 1/2, which disagrees with the evidence, experimental and 
theoretical, that is about 1/3. 

As we showed above, the approach that we propose preserves the mean-field 
critical exponents; all the critical mean-field exponents are suitable. Here we do not 
renormalize the critical exponent from the renormalization-group viewpoint, but 
we take into account the explicit fact that the mean-field critical temperature T c 
is in reality a characteristic scale of temperature linked to the thermal fluctuations 
rather than to the transition temperature. 



4.1.1 Renormalized Weiss theory 

Ferromagnetism and the Curie temperature were explained by Weiss in terms 
of a huge internal " molecular field" proportional to the magnetization. The theory 
is applicable both to localized and delocalized electrons. No such magnetic field 
really exists, but it is a useful way of approximating the effect of the interatomic 
Coulomb interaction in quantum mechanics. When the distance between magnetic 
moments is small, the Pauli exclusion principle, which states that two identical 
fermions may not have the same quantum states, results in interaction between 
magnetic moments. Heisenberg introduced a model to describe this exchange 
interaction on microscopic scale. The Heisenberg exchange Hamiltonian may be 
written in the form 

H exc h = —2 2J JijSi ■ Sj (56) 

i<j 

the summation extends over all magnetic moment pairs in the crystal lattice. 
For positive values of the exchange constant Jij one finds parallel alignment else 
antiparallel. Ferromagnetism is observed for positive exchange interactions below 
a critical temperature. 

The exchange interaction can be regarded as effective field acting on the mo- 
ments. This field is produced by the surrounding magnetic moments and called 
here "renormalized molecular field". As the size of the surrounding moments is 
proportional to the magnetization, the renormalized molecular field is written 
as 

H* m = N* w M } N ^ = N W (1-^) (57) 

with Nw the renormalized Weiss-field coefficient. It is necessary to note that Nw 
was already introduced in the early 20 th century long before the development of 
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quantum physics jUJ H5] . Within the framework of the Weiss theory, one postu- 
lates the existence of a mean field by a phenomenologic approach. The structure 
of parameter iVjy, in such a theory, is not specified. We propose within the frame- 
work of this theory to make a microscopic justification of it; we place ourselves 
within the framework of the "localized magnetism" where the elementary mag- 
netic components are localized on each site of the crystal lattice. The anharmonic 
fluctuations tend to generate mechanisms which are opposed to the local order 
leading to a frustrated internal magnetism. The Weiss coefficient now takes into 
account both the dimension of the system, its size and critical fluctuation and the 
total magnetic field experienced by a magnetic material is thus the sum of the 
externally applied field H and the internal field, 

H = H + H* m , (58) 

and Eq. (54) needs to be rewritten as 

C 1 

X m = — 7^n, where 7 = 1 and T* = T C tt dc . (59) 

1 1 — i c | a 

T* is the renormalized critical temperature and T c would be the critical tempera- 
ture in absence of the correction term. 

For most materials, the phase transition from the paramagnetic state to the fer- 
romagnetic state is found to be of second order. This means that the temperature 
dependence of the first derivative of the free energy is continuous and the second 
derivative of the free energy is discontinuous. Within a "renormalized molecular 
field" model we introduce an additional dimension ^-dependent parameter 

T* c =f{v d )-T c . (60) 

The size dependence of T* enters in the amplitude prefactor f{y d ) defined by 

f(v d ) = l — u d , where v d = (61) 

4.1.2 The localized magnetism 

The (localized) ferromagnetism is characterized by the short-range exchange 
interactions [46J. In the Ising-like model, it is shown for instance that £o ~ d 
(do is the lattice periodicity ), which is in the order of some few A [57]. Indeed, 
the critical width Ate ~ 10 _1 is appreciable and implies a degree of resolution 
experimentally accessible and lets us predict strong corrections to the MFT and 
SGA [9]. Indeed, the discrepancy between the mean- field results and experiment 
though the disagreement between the critical exponents for different dimension- 
alities suggests that the mean-field results are too universal and suggests some 
essential dependence of the critical point on d [38J . To understand the mechanism 
by which £ sets the characteristic length scale of fluctuations we have to consider 
the scattering amplitude. 
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4.2 Ferroelectric thin films 



Finite-size effects in dipolar (magnetic or ferroelectric) ultrathin films is 
a topic of growing interest for both technological and fundamental reasons, (see, 
e.g., references [39-47]). One particularly important issue in these low- dimensional 
systems is the dependency of their Curie temperature (that is, the highest tem- 
perature at which a spontaneous polarization exists) on the film's thickness. The 
technological relevance of knowing such dependency partly stems from the possi- 
bility of designing improved devices since many properties fundamentally depend 
on the Curie temperature. Finite-size scaling theory predicts that this critical 
temperature shifts to lower temperatures than that of the bulk when one or more 
of the material's dimension(s) is reduced to an atomic size [39]. The critical tem- 
perature T C (L) at which the finite system begins to have nonzero average order 
parameter is given by [56] 

T C (L) =Tc-^- (62) 

The existence of a minimum reduced length L m ; n [57] 158] is an indication of the 
first shift of critical point for finite-size system with respect to the bulk one. If 
the length L of the system goes to infinity (thus the system becomes unlimited), 
T C (L) in the last equation returns to the bulk value T c . But for small L the first 
shift of critical temperature will be important. 

Crucial characteristics of an important class of materials are currently unknown. 
In particular, one may wonder if relation among Eq. (62) better describes the 
transition temperatures of ferroelectric thin films. In case that Eq. (62) holds 
in ferroelectric films (for any thickness or "only" above some critical thickness), 
one may also wonder if the mean-field critical temperature, T c , is consistent with 
a given universality class since there has been an intense debate for many years 
on whether ferroelectrics belong to some universality classes or not. This paucity 
of knowledge may arise from the fact that many effects can affect the intrinsic 
T c (L)-versus-L curve in ferroelectric thin films. Examples of such effects are the 
increasing importance of depolarizing fields when decreasing the film's thickness 
[5S] or the increasing release of the strain arising from the substrate for thicker 
films. [BQ]. 

Taking into account Eq. (34) recognizing that the system size L "scales" with 
the correlation length of the bulk system specially in the ordered phase, Eq. (62) 
is rewritten as 2 

T c *(i )=T c *[l- ]. (63) 

If the length l of the system goes to infinity (thus the system becomes unlimited), 
T*(/ ) in the last equation returns to the real bulk value T*. But for small Iq the 
second shift (due to finite-size effects) of critical temperature will be important. 
Specially, if Iq — > ir, T*(/ ) —> 0, so that the limited system will stay in the 
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disordered phase until the temperature becomes absolute zero (It should be kept 
in mind that T* ^ for 2D- and 3D-system [I]. Fluctuations in two dimensions 
are not expected to destroy the ordering as in the case of an isotropic Heisenberg 
ferromagnet because of the anisotropy of the coupling.). Thus there exists a lower 
bound for the size of the system above which there may exist phase transition. 

A deviation of the T* (L)-versus-L was experimentally found in Refs [ST| 162] for 
magnetic films below a critical thickness (note that such deviation carries impor- 
tant information, such as, e.g., the possible predominant role of fluctuations and 
surface on the transition temperature of ultrathin films). 

The traditional wisdom for ferroelectric thin films has been that for film thick- 
nesses smaller than 100 nm the depolarization field will destroy any switchable 
polarization making small particles or thin films non- ferroelectric [5*7] [5*5] . The 
general prediction was that the transition temperature T c will decrease with de- 
creasing size and ferroelectricity will vanish below a minimum critical thickness. 

Recent thin film experiments [63] have however shown that switchable ferro- 
electric films can be made down to 0.9 nm for a crystalline Langmuir-Blodgett 
deposited random copolymer of vinylidene fluoride with trifluoroethylene, (PVDF- 
TrF70:30) on graphite. The minimum thickness is just two mono-layers. The crit- 
ical size for small spherical lead zirconate-titanate (PZT) particles was calculated 
to be 25 A[M]. Similarly it has been reported that switchable ferroelectric films 
can be made of PZT down to 3 or 4 nm [65] in disagreement with earlier theories 
but in agreement with recent theoretical calculations. This shows that finite-size 
effects impose no practical limitation on thin film ferroelectric memory capacitors 
though for some designs tunneling currents may become too large. The lateral 
width of the memory element as well does not represent a limitation at present as 
no change in the coercive field has been observed when the lateral size of the PZT 
cell was decreased from 1 /x to 0.1 //. The voltage necessary for polarization rever- 
sal in such a 100 nm x 100 nm x 100 nm cell is in the technologically accessible 
range of 5 V [55]. The corresponding hysteresis curve has been measured via the 
atomic force microscope (AFM) in the piezoelectric mode [67] . 

Taking into account Eq. (16), Eq. (63) can be put in the following form 

With regard to this last equation, a competition is to be envisaged between vari- 
ous types of fluctuations (thermal fluctuations and fluctuations due to finite-size 
effects) and its consideration in the self-consistent theory can show that finite-size 
effects impose no practical limitation on thin film ferroelectric, in agreement with 
recent theoretical calculations and certain effects observed by now [66J. 

It should be mentioned that finite-size effects such as depressions of T c , and 
reductions in Pg (the spontaneous alignment of dipoles in the ferroelectrics) , have 
indeed been observed in some nano-crystals as small as 250 A in diameter [68J . It is 
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however known that at such sizes this is not an electrostatic phenomenon but seems 
to be due to surface strains or to inhomogeneity effects [HS] • One should also stress 
that the local polarization P as a function of the depth z may actually increase as 
the surface is approached as in the case of PZT [70] leading to an increase in T c , as 
the thickness d decreases. Alternatively it may decrease resulting in a depression 
of T c , with decreasing d. Both effects have been observed by now. This can be 
described by the fact that the extrapolation length can have either sign leading to 
an increase of a decrease of P at the surface [66J. 



4.3 Superconductors 

In renormalized GLW theory the macroscopic wave function of the super- 
conducting state 

■0(r) = i/; (r)exp[— i<p(r)] (65) 

serves as the order parameter with the amplitude squared l^ol 2 = n s being the den- 
sity of the superconducting particles. We use here the macroscopic wave function 
that is a characteristics of the superffuid state of helium and of superconductivity. 
Using the variation method, the two GLW renormalized equations 

— f-iftV + 2eA) 2 ip + a*ip + b*M 2 i) = 0. (66) 
2m 

Js = _(^V^-^V^)-— |^| 2 A. (67) 
m m 

are found. In the first term of Eq. (66), magnetic-field effects are included by 
making the usual replacement V — > V + 2eA. For the charge of the Cooper pairs 
here we have written -2e, although in the original formulation of the theory the 
electronic charge -e was used. 

If we consider a specimen with dimensions much greater than the penetration 
depth, the magnetic field vanishes inside the sample, and ip = const. Therefore, 
only the last two terms of Eq. (66) are relevant, leading to \ip\ 2 = —a* /2b*. Of 
course, this result is identical to Eq. (44). Inserting this expression in Eq. (67), 
we obtain for the current density 

2e 2 \a*\ A , . 

m o* 

Obviously, this expression for the supercurrent is identical to the second London 
equation. However the statistically appearing (thermal and voluminal) fluctuations 
result in an additional current density (or additional electrical conductivity). 

Taking into account Eq. (12), the last equation can roughly be written in the 
form of two contributions: 

js = jGL 4" Jfluct • (69) 
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The first one, jgl, just reproduces the London expression without fluctuations. 
The second term, i.e., the fluctuation part of the supercurrent 

2e 2 

jfluct ^ ^ Q dA (70) 

has a more sophisticated nature. In order to carry out the jfl uct term the an- 
harmonic contributions in the GLW functional, originating from the fourth order 
term, have to be taken into account. Details will be discussed shortly. Conse- 
quently, the GLW renormalized expressions for the penetration depth A* and the 
coherence length £glw are > respectively, given by 



A "^/^R (71) 

^ GLW ^2m\a*\ 

Both characteristic length scales, A* and £glw> have the same dependence on a*. 
Since their ratio ^glw — ^*/£glw * s a function of b* only, it is temperature, 
dimension and size dependent, given by 



, m 2 b* / ©d 

Practically, the mean-field parameter k is the so-called GL parameter that allows 
a distinction between type I and type II superconductors. This could not be 
true in high-temperature superconductors if we must take into account both the 
fluctuations and finite-size effects. The penetration depth, the coherence length, 
and the critical fields are intimately connected. 

As we saw, T c is a characteristic scale of temperature related to thermal fluctua- 
tions and finite-size effects rather than the transition temperature, that also means 
that thermal fluctuations of the equilibrium state also exist above the real critical 
temperature T*. In the normal conducting state the deviation from equilibrium 
can lead to the transient appearance of the superconducting state within certain 
regions, i. e., to the formation of "puddles" of Cooper pairs. These deviations 
from equilibrium are not stable, and they will disappear more or less quickly. The 
statistical appearance of Cooper pairs will become more and more rare the higher 
the temperature, since with increasing temperature the normal conducting state 
becomes more and more stable compared to the superconducting state. There- 
fore, with increasing temperature, larger and larger deviations from equilibrium 
are needed to generate the superconducting state. If we note further that the 
puddles of Cooper pairs represent perfectly conducting regions, we understand im- 
mediately that already above T* due to the fluctuations in the normal conducting 
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state the statistically appearing puddles of Cooper pairs result in an additional 
electrical conductivity, which must strongly increase on approaching T*. 

This influence of the thermal fluctuations can be clearly detected for a number 
of superconductors. In |71j, for example, it shows the transition curve of a bis- 
muth film near the mean-field critical temperature T c . One sees there clearly that 
the full normal resistance is reached only at temperatures considerably above T c 
and the electrical conductance is plotted instead of the resistance, the additional 
conductance a* of the Cooper pair puddles, statistically appearing and vanishing 
again, is particularly clearly visible. 

The additional conductance due to the Cooper pairs (that is well confirmed by 
experiment and that is due to the fluctuations.) can be calculated from the existing 
theories of superconductivity in combination with the theory of fluctuations [T2| 

Eg. ' 



4.3.1 Paraconductivity for T > T* 

The fluctuations of the order parameter for T > T* can also contribute to the 
electrical conductivity, and the pretransitional rise in the conductivity is referred to 
as paraconductivity. The Fourier transform of the nonlocal electrical conductivity 
can be calculated with the aid of the Kubo formula 

i r+oo 

a*(k) = — — / dt(j k (t)U0)) (74) 



2k B T J_ 

where the bracket implies a combined quantum mechanical and statistical average. 
Using the plane wave decomposition of ip(r) given in Eq. (3) and the definition 
of the GLW current operator given in Eq. (67) without magnetic-field effets, we 
have ^_ 

•?'( r ) = — ^2^ q + l')^^' ex P[^ - <J) ■ r] (75) 
Fourier transforming Eq. (75) in exp(iq ■ r) yields 

j(k) = -Y / (2q + k)r q ^ q+ k (76) 
q 

Substituting Eq. (76) into Eq. (74) and limiting ourselves to the k = 0, we obtain 

A 2fc2 r+oo 

= J ^ dt^ q q'(\Ut)\ 2 \M0)\ 2 )- (77) 

We assume terms with q ^ q' are statistically independent and average to zero; 
thus 

A P 2 h 2 r+°° ^ 

= i^f-i J m dtj2m\(r g m q m\ 2 (78) 
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if we assume an exponential decay, the correlation function entering Eq. (78) may 
be written 

«(W*(0)> = -rr 2 * 2 eM~t/r q ), (79) 

h V + S GLW 

which reduces to the Ornstein-Zernicke form Eq. (32) for t = in the case of 
superconductors. The relaxation time, r q will be calculated from the Landau- 
Khalatnikov model. Note that while we had an equation of motion involving the 
phase of the order parameter, we did not introduce an equation of motion for the 
magnitude, |^|. The time scale for achieving phase variations is necessarily slow 
(hydrodynamic) whereas |^| may change rapidly (by the conversion of supercon- 
ductor into normal conductor and vice versa). An exception is near the second 
order phase transition where the free energy becomes "soft" with respect to a vari- 
ation of the magnitude of the order parameter and the relaxational dynamics slows 
down. This intuitive idea was quantified by Landau and Khalatnikov (1954) in 
their discussion of relaxation phenomena in superfluid A He near the lambda point; 
they made the ansatz [73] 

where T is a rate constant (T > 0). The solution of Eq. (80) is given by 

2m 1 

i/) q (t)=ip q (0)exp(-t/T q ), where r q = ^ - ; . (81) 



™ <i 2 + e GLW 

Thermodynamic systems have characteristic microscopic (rapid) relaxation times 
and the relaxation of the vast majority of their internal degrees of freedom proceeds 
on these time scales. There are two important exceptions: (i) modes involving 
degrees of freedom for which conservation laws exist; and (ii) additional modes 
involving a broken symmetry of the system. Additional equations of motion exist 
when the system spontaneously breaks some symmetry. Taking into account Eq. 
(81) and inserting Eq. (79) into the k = limit of Eq. (78), carrying out the 
integration over time, and writing qf = ^q 2 , we obtain 

The integration over negative time is accomplished by replacing exp(— t/r q ) with 
exp(-|t|/r ff ). 

Inserting the expression for T obtained within the framework of time-dependent 
GLW renormalized theory in the vicinity of the renormalized critical temperature 
[Z21E1] 

a T . (83) 
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and taking into account the dependence of the correlation length £q LW on the 
temperature using Eqs. (12), (25) and (72), the singular component expression of 
the additional conductance is obtained as 



h V a T c J V a T ( 

h 



where = 2l d ^d/2) 2 Jo + °° (i+x^) 3 * s a cons t an t integral which is dimension depen- 
dent and £olvk(0) = ^ T ■ In low dimensions d < 4, we set x = q£ to render 
the fluctuation integral, which is UV-finite, dimensionless. 

In Eq. (84) we have indicated the additional conductance for d-dimensional 
samples. 

• In the case of three-dimensional superconductors, the thickness, the width, 
and the length of the sample are large compared to £gw(0). 

^ (r) " 32^ CH ,,(0) i 1 '^)i f + ^) ' <85) 

• For a thin film with length L small compared to £gw(0), we m ay neglect the 
fluctuations along the film normal and the system becomes effectively 2D. 



• In the case of one-dimensional superconductors, the length of the sample is large 
compared to £gw(0), the thickness t and the width g are small compared to 
Z,glw(0) 

tide is the the solution of the self-consistent Eq. (30) for one-, two- and three- 
dimensional sample. Qualitatively, we can easily understand that the sample di- 
mensions must influence the magnitude of the fluctuations, since the Cooper pair 
density can vary only on a length scale of about £,glw- More rapid spatial varia- 
tions require relatively high energies and, hence, practically do not appear. Within 
a sample that is large in all three spatial directions, the Cooper pair density can 
vary spatially in all directions. All these possible configurations must be taken into 
account in the calculation of the additional conductance. For a two-dimensional 
sample, along the shortest extension the Cooper pair density is always constant 
spatially. Hence, averaging over all possible spatial configurations of the Cooper 
pair density along this direction is not necessary. For a one-dimensional sam- 
ple, averaging is unnecessary along both directions in which the sample is small 
compared to ^glw- We see that the statistics is restricted because of the sample 
geometry. This results in various expressions for the additional conductance. 
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Experience shows that the transition curves of three-dimensional samples, say, 
of wires with a diameter large compared to £q LW , are very sharp, i. e., the effects 
we have just discussed cannot be observed. The reason is not the absence of 
fluctuations, but rather the comparatively high residual conductance of the three- 
dimensional sample. The quantity /(f^c) = — SoT~) °^ ^Q s - (85-87) is an 
indication of the decreased additional conductance as the fluctuation effects of the 
system increase. 

So far we have only discussed how the fluctuations affect the electrical conduc- 
tance. However, if puddles of Cooper pairs appear statistically above T*, this must 
be noticed also in other properties. We know that below T* a superconductor ex- 
pels small magnetic fields out of its interior, i. e., it turns into an ideal diamagnet 
[75] . We expect that, similar to the effect of the fluctuations on the conductance, 
some part of this diamagnetic property also appears above T*. The puddles of 
Cooper pairs should result in a characteristic temperature dependence of the dia- 
magnetic behavior of the superconductor above T*. Only a few hundredths of a 
degree away from T* the additional diamagnetism is already very small and corre- 
sponds to the expulsion of just a few flux quanta. However, it has been possible to 
detect this effect clearly [73] by utilizing a superconducting quantum interferome- 
ter. The fluctuations should also lead to an increase of the specific heat C already 
above T*. This effect could also be experimentally demonstrated [73"| 175]. 

4.3.2 Conventional and unconventional superconductors 

The considerable success of the mean-field theory in conventional supercon- 
ductors originates from the low value (Ate = Qdc/aoT c ~ 10~ 14 ) of the critical 
Ginzburg width [1, 18, 45]. Accordingly, Eq. (16) obviously leads to T* ~ T c , 
where T c is the mean field (BCS) value of the transition temperature. This implies 
that it is impossible to detect with actual experimental precision the deviations 
from the mean-field theory. Indeed, one observes that the specific heat does not 
present any divergence, but rather a compatible jump with the predictions of 
Eq. (3) in [I]. Practically, critical phenomena should not be observable in usual 
conventional superconductors, which could not be true in high-temperature super- 
conductors and particularly high-T c superconductors [3"t 177]. 

In compounds of the BiSrCaCuO, LaSrCaCuO and TISrCaCuO type, estima- 
tions around £ ~ 1 nm were reported, the paired electrons of the Cooper pair (the 
underlying unit for superconductors) are forced apart by their Coulomb repulsion, 
resulting in a relatively large separation of £o ~ 10 3 A [31 H] • These values have 
to approximately equal the size of the units that undergo ordering at the phase 
transition. This can question the validity of the mean-field theory for supercon- 
ductors. But, we think that such a distance is related to the existence of large 
fluctuations and finite-size effects present in those compounds. The increase of 
the critical temperature in non-conventional superconductors also means that the 
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bond energy of the Cooper pair increases at the same rate. When the energy bond 
is higher, the bond is mostly confined in space and £o then reduces at the same 
level. It is for instance the case in La>2- x Sr x CuO '4 where Ate can reach the value 
1CT 1 , and in Y Ba2Cu^,0-j- y compounds where Ate can reach the values typically 
in the range 1CT 2 to 10 _1 [4j. These values seem to be accessible in experiment 
and suggest principal deviations compared to the mean-field theory observable. 
Also, the newer ceramic high-temperature superconductors have a much smaller 
coherence length of £0, which can indeed help to show some effects of fluctuations 

IZSl CZSl 

5 Conclusion 

We have described the fundamentals of renormalized 4 theory on the basis of 
a GLW calculations, combined with an efficient no perturbative technique, that 
takes into account both dimension, size and microscopic details of the system, 
and which leads to critical behavior, strongly deviating from the classical MFA far 
from the thermodynamic limit. Within this more rigorous approach the effects of 
fluctuations are examined in more rigorous detail, and we are able to establish the 
insufficiencies of the MFT, and also estimate the width of the critical region where 
corrections to MFT are important. We find a consistent interpretation relating the 
dimension, the size of the system and the spatial fluctuations which can give an 
interpretation to the degeneracy of energy levels. We have also calculated asymp- 
totic expressions of thermodynamic observables as a function of temperature. 

This theory does not present any renormalization based method in the sense of 
a renormalization group approach but rather a SCM improving on LT. The SCM 
is a strategy for dealing with problems involving many length scales. The strategy 
is to tackle the problem in steps, one step for each length scale. In the case of 
critical phenomena, the problem, technically, is to carry out statistical averages 
over thermal fluctuations on all size scales. The SCM is to integrate out the fluctu- 
ations in sequence, starting with fluctuations on an atomic scale and then moving 
to successively larger scales until fluctuations on all scales have been averaged out. 
The integration of this deviation in the self-consistency on the entire spectrum of 
the lattice vibration, conducted to a very improved self consistent energy. This 
improvement of the self consistent problem by the processed harmonic and anhar- 
monic fluctuations goes for example to reflect itself on the thermodynamic and 
electronic parameters, resulting notably in a substantial reduction of the critical 
temperature. 

The resulting self consistent problem of this approach has been solved ana- 
lytically, what permitted us to extract an effective theory, with notably a mean 
field critical temperature renormalized by fluctuations. An interesting point with 
this renormalized critical temperature is its dependence on the quantity ^ which 
according to the text, characterizes the importance of the fluctuations and is a 
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function of both dimension and size of the system. 

We have revealed the thickness dependence of the Curie temperature for various 
types of thin and ultrathin films. In agreement with recent theoretical calculations 
it is found that the relation of Eq. (64) reproduces rather well such dependence for 
thickness down to Iq, in contrast with the usual finite-size scaling law that breaks 
down below L and that does not take into account both the fluctuations and the 
microscopic details of the system. We are confident that the present work provides 
a deeper knowledge of nanoscience, phase transitions, and critical behaviors in 
dipolar systems. Furthermore, the parameters characterizing this approach depend 
strongly on the dimension and the size of the sample, but in different ways than 
the ones obtained in the MFA. Thus, this SCM seems to be appropriate for model 
calculations in a number of critical systems. For this reason, analysis has also been 
performed on other systems such as non-conventional superconductors, localized 
magnetism and can be extended to all system, whose transitions indeed apparently 
belong to the same universality class, in order to lend a significant level of support 
to this theoretical approach. 

By comparing the thermal energy and the energy resulting to structural insta- 
bilities, we introduced the report which could allow to estimate the quantum 
fluctuation importance. For large values of u d , only a quantum description can give 
account of the dynamics of the system and we must use another approach that 
is the second step in this theory and that is to write down a general field theory 
for the order parameter, consistent with all symmetries of the underlying model. 
As we are dealing with a quantum transition, the field theory has to extend over 
space time, with the temporal fluctuations representing the sum over histories in 
the Feynman path integral approach. 

Moreover, the SCM appears very suitable for providing a quantitative descrip- 
tion in the crossover region as its application makes it possible to address the 
importance of fluctuations and thence to elucidate the basic physics of phase tran- 
sitions and critical phenomana. Finally, we also expect that our method will reveal 
universal features of critical systems with finite size, which needs numerical studies 
and pronounced experimental confirmation. 
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